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ABSTRACT 

We analyze measured proton and electron temperatures in the high-speed solar wind in order to calculate the 
separate rates of heat deposition for protons and electrons. When comparing with other regions of the helio- 
sphere, the fast solar wind has the lowest density and the least frequent Coulomb collisions. This makes the 
fast wind an optimal testing ground for studies of collisionless kinetic processes associated with the dissipation 
of plasma turbulence. Data from the Helios and Ulysses plasma instruments were collected to determine mean 
radial trends in the temperatures and the electron heat conduction flux between 0.29 and 5.4 AU. The derived 
heating rates apply specifically for these mean plasma properties and not for the full range of measured values 
around the mean. We found that the protons receive about 60% of the total plasma heating in the inner helio- 
sphere, and that this fraction increases to approximately 80% by the orbit of Jupiter. A major factor affecting 
the uncertainty in this fraction is the uncertainty in the measured radial gradient of the electron heat conduction 
flux. The empirically derived partitioning of heat between protons and electrons is in rough agreement with 
theoretical predictions from a model of linear Vlasov wave damping. For a modeled power spectrum consist- 
ing only of Alfvenic fluctuations, the best agreement was found for a distribution of wavenumber vectors that 
evolves toward isotropy as distance increases. 

Subject headings: hydrodynamics — MHD — plasmas — solar wind — turbulence — waves 



1. INTRODUCTION 

The supersonic solar wind is accelerated away from the Sun 
by some combination of physical processes including gradi- 
ents in gas pressure (from the hot, 10 6 K corona) and wave 
pressure, as well as possible collisionless wave-particle inter- 
actions. There is gradual heat input into the solar wind plasma 
that begins in the corona and extends far out into interplane- 
tary space (see reviews by Parker 1963; Leer et al. 1982; Tu 
& Marsch 1995; Goldstein et al. 1995; Marsch 1999; Holl- 
weg & Isenberg 2002; Cranmer 2002; Matthaeus et al. 2003). 
One likely source of this extended heating is the dissipation of 
magnetohydrodynamic (MHD) turbulence. In order to better 
understand the combined problem of coronal heating, solar 
wind acceleration, and the large-scale evolution of the tur- 
bulent heliospheric plasma, we need to know how energy is 
transferred between the MHD fluctuations and the particles. 

In the mainly collisionless solar wind, the various parti- 
cle species (i.e., protons, electrons, and heavy ions) are not 
in thermal equilibrium with one another. The particles ex- 
hibit a range of different outflow speeds, temperatures, and 
velocity distribution anisotropies, and these differences are 
most pronounced in low-density regions with the least fre- 
quent Coulomb collisions (e.g., Neugebauer 1982; Kohl et al. 
1997, 1998, 2006; Kasper et al. 2008). These differences can 
be used to probe the kinetic physical processes that are re- 
sponsible for depositing energy into the plasma. 

There have been a number of studies where the measured 
plasma properties in the solar wind were used to derive the 
corresponding rates of energy input from processes such as 
MHD turbulence (e.g., Coleman 1968; Tu 1988; Freeman 
1988; Verma et al. 1995; Matthaeus et al. 1999b; Smith et 
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al. 2001; Vasquez et al. 2007; MacBride et al. 2008; Marino 
et al. 2008; Breech et al. 2008). Much of this work, however, 
dealt only with the energy budget of the protons and not with 
any other ions or the electrons. Although the solar wind's 
mass density and momentum flux are dominated by the pro- 
tons, the electrons carry approximately half of the thermal en- 
ergy of the plasma and should not be neglected in a complete 
treatment. There have been other investigations into the elec- 
tron energy balance in the solar wind (e.g., Scudder & Olbert 
1979; Phillips & Gosling 1990; Pilipp et al. 1990). There also 
have been studies of turbulent dissipation that compare a total 
heating rate to the presumed proton contribution — and thus 
estimate the electron heating rate as the remainder (Leamon 
et al. 1999; Stawarz et al. 2009). However, there has been sur- 
prisingly little combined analysis of the proton and electron 
heating rates that treat the two plasma components on equal 
footing. 

In this paper we compute new estimates of proton and elec- 
tron heating rates in the fast solar wind. These rates are 
derived from plasma properties measured by the Helios and 
Ulysses spacecraft. In § 2 we describe these measurements in 
detail, and in § 3 we show how the heating rates can be com- 
puted from the separate equations of energy conservation for 
protons and electrons. We find, not surprisingly, that an ac- 
curate determination of the electron heating rate depends cru- 
cially on the measured electron heat conduction flux. In § 4 
we compare the empirically derived heating rates with heat- 
ing rates obtained from linear Vlasov theory, in which several 
simple assumptions about the symmetry of an (Alfvenic) fluc- 
tuation spectrum were made. Finally, § 5 contains a summary 
of the major results of this paper and a discussion of the im- 
plications these results may have on our wider understanding 
of heliospheric plasma physics. 

This work is being presented in tandem with an indepen- 
dent investigation by Breech et al. (2009). In both papers, we 
studied the same problem of proton-electron heat partition- 
ing in the fast wind, but it has been approached from differ- 
ent and complementary vantage points. This paper attempts 
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to "invert" the in situ measurements in order to determine the 
partition fraction between proton and electron heating. On the 
other hand, Breech et al. (2009) showed how a sophisticated 
model of MHD turbulent heating — with an assumed value for 
the partition fraction — can also be consistent with the empir- 
ical data. The results of these two studies are consistent with 
one another. 

2. IN SITU PARTICLE DATA 

We examine proton and electron plasma properties for the 
high-speed solar wind between 0.29 and 5.4 AU. Measure- 
ments made at larger distances (e.g., from the Voyager probes) 
are excluded for two reasons: (1) these data have been taken 
mainly in the ecliptic plane, and thus are dominated by slow- 
speed wind, and (2) the internal energy budget of particle ve- 
locity distributions at distances r > 10 AU is increasingly af- 
fected by pickup ions, which are ignored here. We focus on 
fast wind streams because these appear to be the sites of the 
most ambient and "quiescent" solar wind plasma, and because 
their low densities highlight them as the least complicated by 
inter-species collisional coupling. 

Figure 1 displays the data we use for the proton temperature 
T p , electron temperature T e , and electron parallel heat conduc- 
tion flux qn e . Operationally, the high-speed wind was defined 
as being all streams faster than 600 km s" 1 . This selection cri- 
terion was used for both the Helios and Ulysses data sets. A 
speed of 600 km s" 1 is slightly more restrictive than the more 
standard value of 500 km s" 1 that has been applied by others 
to define the fast wind (e.g., Dasso et al. 2005; MacBride et 
al. 2008). Taking a narrower range of speeds results in less 
contamination from plasma parcels with qualitatively differ- 
ent properties. 

All measurements for heliocentric distances r less than 1 
AU came from published data from the plasma instruments 
on the two Helios spacecraft (Schwenn & Marsch 1990). 
Proton temperatures for the fast wind were obtained from 
Marsch et al. (1982) and electron plasma properties (T e and 
q\\ te ) were taken from Pilipp et al. (1990). The tabulated data 
points, which spanned the near-solar- minimum years 1974- 
1976, corresponded to representative intervals that were se- 
lected for broad coverage in distance and wind speed. The 
actual velocity distributions for both protons and electrons 
were non-Maxwellian, and the temperatures used here are 
mean isotropic (drifting Maxwellian) equivalents, i.e., T = 
(7||+27i)/3. 

Data for r > 1 AU came from the "Solar Wind Observations 
Over the Poles of the Sun" (SWOOPS) instrument on Ulysses 
(Bame et al. 1992; Goldstein et al. 1996). The actual temper- 
ature measurements were extracted from the European Space 
Agency's online archive. 4 We collected all SWOOPS data for 
the decade of time starting at the beginning of the Ulysses mis- 
sion (1990 November) and extending to 2000 December. This 
fully encompassed the solar minimum in 1996-1997, during 
which there were found to be large swaths of the heliospheric 
volume containing quiescent fast solar wind. The highest time 
resolution data available for each interval was used to better 
pinpoint regions of high-speed wind (using the proton outflow 
speed u p ). For the proton temperature data, this corresponded 
to a 4 or 8 minute cadence, depending on the mode of op- 
eration. For the electron temperatures, the highest cadence 
tended to vary between 7 and 35 minutes. 

4 The ESA Ulysses D ata Syste m (UPS) site for SWOOPS data is: 
http://helio.esa.int/ulysses/archive/swoops.html 
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FIG. 1 . — In situ measurements for the fast solar wind: (a) plasma temper- 
atures from Helios (filled diamonds) and Ulysses (small points), with protons 
in red and electrons in blue, (b) electron heat conduction flux from Helios 
(filled diamonds) and Ulysses (triangles). Also shown are least-squares fits 
(solid lines), the classical Spitzer-Harm heat flux (dashed line), and the Holl- 
weg (1974, 1976) collisionless heat flux with a e = 1.05 (dotted line). 



The Ulysses electron heat conduction data shown in Figure 
lb were taken from Figure 7 of Scime et al. (1994). These 
measurements were taken in the ecliptic plane in 1990-1992 
during the initial Ulysses cruise phase from 1 to 5 AU. These 
data were not provided in a format where the high-speed wind 
selection criterion (u p > 600 km s -1 ) could be applied. Thus, 
the points in Figure lb at r > 1 AU apply to all wind speed 
intervals. However, Scime et al. (1999) found that there was 
no significant variation in qu e as Ulysses passed between fast 
and slow wind streams during its first high-latitude scan in 
1994-1995 (see also Salem et al. 2003). 

There is a drastic difference in the number of data points 
shown in Figure la between the Helios and Ulysses mea- 
surements. For Helios, the high-speed velocity criterion re- 
sulted in only 21 data points for T p and 18 data points for T e , 
For Ulysses, however, there were 211593 T p data points and 
61865 T e data points for the fast wind. This is actually only 
a subset of the total amount of data present in the archive, 
since we used only proton data validated by a quality control 
flag that indicated nothing detectably wrong with the mea- 
surement. In order to keep Figure 1 from being unreason- 
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ably large in size, we plotted every 20th individual SWOOPS 
temperature measurement. For the collection of Ulysses data 
(chosen by the criterion of u p > 600 km s" 1 ), the mean wind 
speed was found to be 744 km s" 1 , with a standard deviation 
about this mean of 47 km s . If the data selection criterion 
is reduced to u p > 500 km s" 1 , the mean speed and standard 
deviation change to 705 km s" 1 and 92 km s" 1 , respectively. 

There are a number of potential difficulties and systematic 
uncertainties in measuring T p with the Ulysses plasma instru- 
ment. Some of these issues are more problematic for the slow 
wind, since these regions tend to have lower values of T p and 
thus narrower distribution functions in velocity space. These 
measurements may suffer from discretization effects caused 
by the finite number of energy and angle channels of the 
SWOOPS instrument. However, like many solar wind instru- 
ments, SWOOPS has a logarithmic distribution in the spac- 
ing of its energy channels. This could result in an additional 
source of error for the fast wind, since the velocity distribu- 
tions are "centered" at higher speeds where the relative energy 
resolution is coarser and the protons are thus spread across 
fewer energy windows. 

In order to convey a sense of the measurement uncertainties, 
the SWOOPS data archive reports two independent determi- 
nations of the proton temperature that are claimed to usually 
bracket the true proton temperature from below (r sma u) and 
from above (7]^): 

1 . The lower limit value r sma n is essentially just the ra- 
dial component of the proton temperature, with the 
high-energy "beam" component cut off. Goldstein et 
al. (1996) concluded that the proton velocity distribu- 
tion is roughly isotropic — over the Ulysses distances — 
because there was never a substantial dependence of the 
radial temperature on the local magnetic field direction 
(see also Vasquez et al. 2007). Thus, the closer the dis- 
tribution tends to being isotropic in velocity space, the 
better an assumption T sma a will be. 

2. The upper limit value 71 arge was obtained by integrating 
the measured proton velocity distribution over all of the 
energy channels and angle bins that were statistically 
above a determined noise level. For the lowest tem- 
peratures, this process results in an overestimate of T p 
because the angular responses of the SWOOPS instru- 
ment channels may be broader than the actual velocity 
distribution function. 

In Figure la, we adopt the geometric mean of these two 
bounding values to obtain an intermediate estimate for T p = 

(7smaii7iarge) 1//2 - It should be noted that there seems to be a 
solar cycle trend in the ratio T^ge/T^u. At solar maximum 
(i.e., around 1990 and 2000), this ratio tends to range between 
2 and 3. At solar minimum (i.e., 1996-1997), the ratio tends 
to range between only 1.2 and 1.5. Thus, the more we depend 
on the solar minimum SWOOPS data (see below), the less we 
need to be concerned about the differences between 7i alge and 
Tsmaii (see also Marino et al. 2008). 

For the Ulysses electron data, we used the total tempera- 
tures formed by a weighted sum of the narrow thermal core 
and the broader "halo" and "strahl" components of the veloc- 
ity distributions. The total electron temperature is the most 
consistent quantity to use with the internal energy moment 
equations described in § 3.1, since these equations were de- 
rived by integrating over all velocity space. Figure la shows 



a noticeable bifurcation of T e into hot and cold branches at 
distances around r w 3^4 AU. This appears to be a solar cycle 
effect, since the lower values are more dominant during 1996 
(solar minimum) and the higher values are more favored dur- 
ing the two solar maximum periods (1991 and 2000) and, to 
a lesser extent, the post-maximum phase (1994). The higher 
T e data points at 3-4 AU also exhibited slightly higher elec- 
tron densities n e than the cooler data points. Thus, the lower 
values seem more appropriate to be applied to studies of the 
ambient fast solar wind associated with polar coronal holes at 
solar minimum. The least squares fit curve shown in Figure 
la for T e serendipitously favors these lower values. 

Figure 1 also shows analytic fits for the measured quantities 
between 0.29 and 5.4 AU. These fits are given by 

In (jfip^j = 0.03460-0.4333x+0.08383x 2 (2) 

ln(ik) = -0.7032 -2.1 15jc-0.2545jc 2 (3) 

where x = ln(r/[l AU]) and qo = 0.01 erg cm" 2 s" 1 . To avoid 
the larger number of Ulysses data points overwhelming the 
Helios data points in the least-squares fitting process, the data 
sets from the two spacecraft were weighted equally. These 
fits should not be extended very far beyond the range of he- 
liocentric distances (0.29 < r < 5.4 AU) for which they were 
derived. 

The new fits given above can be analyzed in terms of a com- 
monly assumed power-law dependence of temperature as a 
function of radius; i.e., T oc r~ s . The local value of the ex- 
ponent 5 can be computed at any distance as the logarithmic 
derivative -dlnT /dlnr. Using equations (fl]i and (f2]), the re- 
sulting proton and electron temperature exponents S p and & e 
both decrease monotonically with increasing distance. For 
protons, S p is approximately 0.98 at 0.3 AU, it decreases to 
0.80 at 1 AU, and then to 0.56 at 5 AU. This range of val- 
ues agrees with other exponents reported over these distances 
that range between about 0.75 and 1 (e.g., Eyni & Steinitz 
1981; Lopez & Freeman 1986; Totten et al. 1995; Ebert et 
al. 2009). For electrons, S e is approximately 0.65 at 0.3 AU, 
it decreases to 0.43 at 1 AU, and to 0.15 at 5 AU. These are 
also in agreement with earlier measurements of order 0.2-0.6 
(see Sittler & Scudder 1980; Pilipp et al. 1990; Phillips et al. 
1995; Issautier et al. 1998). Some of these studies did not 
include any explicit criteria for selecting either fast or slow 
solar wind streams, but it is interesting that the same rough 
range of exponents is found. It is possible to use these expo- 
nents to derive an effective polytropic index for the protons 
and electrons — i.e. ,7= l + (<5/2) (Totten et al. 1995) — but the 
analysis in § 3 attempts to do a more thorough study of the 
internal energy balance of the plasma. 

Figure lb also shows how the measured electron heat con- 
duction flux compares to both the classical collisional heat 
flux (Spitzer & Harm 1953) and to Hollweg's (1974, 1976) 
collisionless "free streaming" approximation. The former is 
given by 

q\\,e = -Ke-^ (4) 

where the electron conductivity is 

T 5/2 

K e = (1 .84 x 10" 5 erg cm -1 s _1 KT 7/2 ) . (5) 
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Also, the electron Coulomb logarithm is approximated by 

lnA„ = 23.2 + -lnf-^-V-lnf r 1 * ) ( 6 ) 
2 \1Q 6 KJ 2 V10 6 cm- 3 / V ' 

where n e is the electron number density (see also Cranmer et 
al. 2007). The Spitzer-Harm values are reasonably close to 
the measured data above r f» 1 AU, but they exceed the mea- 
surements at smaller distances. The Hollweg (1974, 1976) 
heat flux is given by 

q\\, e ~ ^a e n e uk B T e , (7) 

where k B is Boltzmann's constant. The dimensionless order- 
unity constant a e is only known approximately, and it is ex- 
pected to depend on the microscopic shape of the electron 
velocity distribution (Hollweg 1974, 1976). Using the above 
expression, though, we solved for a e using each of the data 
points in Figure lb. For simplicity, we used the above fit for 
T e (r), we also assumed that u = 700 km s -1 , and we used the 
empirical analytic model for n e {r) that is described below. The 
mean value of a e for all 38 points was 1.05, with a standard 
deviation about the mean of 0.44 and no clear radial trend. 
(The full set of values ranged between 0.47 and 2.08.) For 
comparison, Figure lb shows a corresponding curve for the 
collisionless heat conduction flux when a e = 1.05. 

3. EMPIRICAL HEATING RATES 
3.1. Internal Energy Conservation Equations 

Our goal is to use measurements to quantify the rates of 
heating and cooling associated with as many physical pro- 
cesses as possible, and then to solve for the net volumet- 
ric rates of heat input that presumably can be attributed to 
MHD turbulence. The primary mechanisms that are taken 
into account include adiabatic energy conservation, electron 
heat conduction, and Coulomb collisions. We also retain the 
assumption that both the proton and electron velocity distri- 
butions are isotropic Maxwell-Boltzmann distributions with a 
shared bulk velocity (u p = u e ). Direct collisions between the 
proton and electron populations are included for complete- 
ness, but they are expected to be extremely weak in the helio- 
sphere. The models computed using classical collision rates 
(e.g., Spitzer & Harm 1953) are virtually identical to those 
without any collisions whatsoever. However, we also study 
the ramifications of more rapid collisions by scaling up the 
classical rates by a range of constant factors. 

It is important to discuss some of the physical processes 
that are neglected in the models presented below. For exam- 
ple, we do not explicitly consider how departures from tem- 
perature isotropy (i.e., 7ji ^ 7j_) affect the energy balance. Ob- 
servationally, these departures are not significantly large over 
most of the distances considered here (see Marsch et al. 1982; 
Pilipp et al. 1990; Matteini et al. 2007). At 1 AU, Vasquez et 
al. (2007) showed that taking account of the measured range 
of anisotropies would result in, at most, about a 10% change 
in the adiabatic terms of the energy conservation equations. 
Still, it has been shown that for some regions in the solar 
wind, the velocity distributions are kept from deforming too 
far away from isotropy by plasma instabilities (e.g., Kasper et 
al. 2002; Hellinger et al. 2006). We defer any consideration of 
the energy transfer between instability-generated waves and 
the particles to future work. 

Another process that is neglected here is proton heat con- 
duction. The classical proton heat conduction coefficient is 



about 25 times smaller than that for electrons (Braginskii 
1965; Sandbsek & Leer 1995). Also, some models based 
on higher-order moment closures of the Boltzmann equation 
have found even lower values than the classical approach 
would suggest (e.g., Olsen & Leer 1996). The proton heat 
conduction flux was measured by Helios between 0.3 and 1 
AU (Marsch et al. 1982). To confirm the expectation that the 
proton heat flux can be neglected, we analyzed these measure- 
ments in terms of the collisionless free-streaming approxima- 
tion discussed above. In other words, for each measurement 
we computed an effective coefficient a p analogous to the elec- 
tron coefficient a e in equation (|7). Using the data tabulated by 
Marsch et al. (1982), we obtained a mean value of a p = 0.045, 
which is a factor of ^23 lower than the mean electron coeffi- 
cient a e = 1.05. Test runs of the energy conservation models 
that included this level of proton heat conduction resulted in 
heating rates that were within 5% of the ones computed with- 
out this effect. 

We are now in a position to describe the conservation laws 
for proton and electron internal energy upon which the present 
results will be based. Assuming a time-steady solar wind that 
has reached a constant asymptotic terminal speed u (see, e.g., 
Arya & Freeman 1991), these equations are 

^n p uk B ^ - uk B T p ^ = Q p + ^n p k B v pe {T e - T p ) (8) 

\n e ukJ^-uk B T e ^- = Q e + \n e k B v ep {T p -T e ) 
2 or or 2 

-^J^H^cos 2 *) (9) 

where the heat input rates are Q p and Q e , the Parker spiral an- 
gle is <I>, and the rates of proton-electron Coulomb collisions 
for the two equations are v pe and v ep (see Barakat & Schunk 
1982; Isenberg 1984; Cranmer et al. 1999). As mentioned 
above, we assumed that the outflow speed u is constant and 
identical for the protons and electrons. We also assumed that 
the electron densities n e and n p vary with distance as r" 2 , and 
that n p is normalized to a value of 2.5 cm" 3 at 1 AU (e.g., 
Goldstein et al. 1996). We utilized a 5% helium abundance 
by number in order to compute n e = \.\n p . Also, the winding 
angle of the spiral interplanetary magnetic field is given in its 
standard form as 

tan$ = ftrsin6»/ M , (10) 

where we used a rotation frequency £1 = 2.7 x 10~ 6 rad s -1 . In 
most of the models shown below, we set the colatitude 9 = 15° 
to model the high-latitude Ulysses measurements. 5 

The Coulomb collision frequencies are balanced such that 
n e u ep « n p v pe . For two isotropic Maxwellian distributions in- 
teracting with one another, 

32 ./,, . e 4 n e n„ „ 

n p y p e — — 7T ' lnA '-^ (H) 

3 m e m p a 3 

where the Coulomb logarithm In A w 27 at 1 AU, m p and m e 
are the proton and electron masses, and e is the magnitude 

5 Although the Helios measurements were made close to the ecliptic plane, 
the Parker spiral effect in the inner heliosphere is not as pronounced as it is at 
r > 1 AU because of the linear dependence on distance. Thus, the resulting 
values of Q e at the lowest values of r are relatively insensitive to the choice 
for 9 (see Fig. 4b). 
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of the proton and electron charge (Spitzer 1962). The mean 
squared interaction speed is given by 

o 2 = M + M. (12) 

m p m e 

Note that, as defined above, the quantity v pe is a factor of 
two larger than the equivalently named rate used by Isen- 
berg (1984) and Cranmer et al. (1999). The present definition 
is more consistent with it being the true rate of temperature 
equilibration as described by Spitzer (1962). After evaluat- 
ing many of the constants and using the approximation that 
m e <C m p , we found 

^^■4xl0-'(^)(^)"' 2 s-'. (13, 

This expression gives rise to a rather large mean free path for 
electron-proton collisions. Depending on whether the mean 
free path is defined in terms of the solar wind speed (L m f p ~ 
// / v pe ) or the faster electron thermal speed (L m f p ~ V,h, e / v pe ), 
this quantity is of order 500-1500 AU at r = 1 AU. This should 
be contrasted with the much smaller mean free path of ~0.5 
AU for electron-electron self-collisions that maintain the ther- 
mal core at 1 AU (see, e.g., Spitzer 1962; Salem et al. 2003). 

In addition to Coulomb collisions, there may be other 
collision-like processes that lead to temperature equilibration 
and isotropization in the heliosphere. For example, collision- 
less wave-particle interactions have been suggested for many 
years as being able to produce these effects (Cuperman & 
Harten 1970; Perkins 1973; Dum 1983; Williams 1995; Kel- 
logg 2000). If the actual proton-electron temperature equili- 
bration rate is faster than expected, then more heat must go 
into the protons for them to maintain the known inequality 
T p > T e in the fast wind. Thus, if the collision rate is anoma- 
lously enhanced, the resulting value of Q p would be larger 
(and Q e would be smaller) than the values computed with 
weak or nonexistent collisions. We include these effects be- 
low by multiplying the collision rates v pe and v ep by an arbi- 
trary constant /. 

The electron heat conduction flux affects the energy bal- 
ance in qualitatively different ways depending on heliocentric 
distance. Using the empirical fit for q\\, e {r) (eq. [3]) and tak- 
ing the divergence as shown in equation (0, we found that 
this term leads to local electron cooling for r < 0.75 AU, be- 
cause the radial slope of q\\ :C is flatter than r~ 2 . In this case, 
cool plasma at larger radii is conducted inward. At larger dis- 
tances, though, this term gives local electron heating, since 
the slope of qu e is steeper than r~ 2 and hot plasma is con- 
ducted outward. Note, however, that if we had used either the 
classical Spitzer-Harm heat flux or the collisionless a e = 1 .05 
approximation, as shown in Figure lb, then qn e would always 
be slightly steeper than r~ 2 and heat would be conducted out- 
wards at all distances. 

3.2. Results for Proton and Electron Heating 

We solved equations © and (0 for the volumetric heating 
rates Q p and Q e over the range of heliocentric distances cov- 
ered by the Helios and Ulysses measurements. An interesting 
aspect of this work is that because the internal energy equa- 
tions are not being solved for the temperatures — but instead 
for the heating rates — we can avoid complicated numerical 
differential equation techniques and use a straightforward al- 
gebraic solution for Q p and Q e . Even the radial derivatives 
can be computed analytically from the fits given in § 2. In 



1 1 — i — i 

o.oi 




0.1 1 10 

r (AU) 



FIG. 2. — Empirically derived heating rates for protons (solid lines) and 
electrons (dashed lines) in the fast solar wind, with multiple curves showing 
results for u = 600, 650, 700, 750, and 800 km s" 1 (from bottom to top for 
each set of curves). Shown for comparison is Q lur b (eq. 1151 ) for five values 
of Aj^ at 1 AU (dotted lines). All heating rates have been multiplied by 
(r/1 AU) 4 . Example error bars are given for the u = 700 km s case (see 
text). 

practice, however, we computed the quantities to be differen- 
tiated on a very fine grid (1000 points between 0.27 and 5.5 
AU) and used a standard centered-difference approximation 
to compute each of the terms in equations ([8]) and ©. 

The internal energy equations were first solved for two 
cases: (1) standard Coulomb collision rates, with / = 1, and 
(2) a completely collisionless heliosphere, with / = 0. These 
two cases gave virtually identical results to one another, with 
relative differences between the heating rates only at the 0.5% 
level. The procedure was also carried out several times for a 
range of assumed (constant) wind speeds u between 600 and 
800 km s -1 . One can see from equation ^ that, in the case of 
negligible Coulomb collisions, our derived proton heating rate 
Q p should be linearly proportional to the wind speed u. The 
electron heating rate Q e also increases as the wind speed in- 
creases, but because of the heat conduction this is not a purely 
linear relationship. 

Figure 2 shows the quantities r A Q p (r) and r Q e (r) as a func- 
tion of radius, rather than the rates themselves, because the 
latter decrease very steeply with distance. It is easier to see 
the subtle relative differences between Q p and Q e when the 
dominant radial variation has been removed. It is clear that the 
computed values for Q p can be well approximated by power 
law scaling relations in both radius and wind speed. The fol- 
lowing fit was found to be valid to within about 6% relative 
accuracy over the full range of distances: 

Q,, w 3.42 x 10" 16 f— — ) 35 | - r ) ergs -1 cm -3 . 

* P VlAU/ 1,700 km s" 1 ; S 

(14) 

These scalings compare favorably to similar calculations by 
Verma et al. (1995) and Vasquez et al. (2007). 

Note that it is not possible to fit the electron heating rate 
Q e {r) with a simple power-law function because of the non- 
local heat conduction. The distances at which Q e decreases 
rapidly (~3 AU) are roughly the same as where the "bifurca- 
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tion" in T e occurs (see Fig. 1). We cannot conclude whether 
this is a coincidence or not because of the relatively sparse 
sampling of the fast wind by Ulysses at these distances, Note 
also that solving equation (0 gives negative values for Q e 
at these heights for the smallest wind speed of 600 km s -1 . 
The relative impact of the electron heat flux is strong at these 
distances, and the negative values would be eliminated if the 
magnitude of qn e was smaller by only about 10%. This is 
certainly within reason, especially since the data used to con- 
strain the fit (Scime et al. 1994) were not specifically taken in 
the fast solar wind. 

Figure 2 also includes representative uncertainty limits, 
shown as error bars only at two distinct heights in order to dis- 
tinguish them clearly from the curves showing the variations 
in u. These uncertainty limits reflect the existence of a range 
of measured values around the mean fitting curves used for T p , 
T e , and qn e . To compute these uncertainty bounds on Q p and 
Q e , we created a set of alternate models in which each of the 
fits was was multiplied by factors 1.3", where n = -1,0, +1. 
The result was a grid of 27 models in which the tempera- 
tures and electron heat fluxes were varied up and down (in 
all combinations) by fiducial ±30% amounts. The extreme 
upper and lower limits on Q p and Q e were found for all of 
these models and are plotted in Figure 2 at the two example 
distances. As can be seen from Figure 1, the actual spread in 
the in situ data often reaches — and sometimes exceeds — this 
fiducial 30% relative variability level. Thus, it is clear that the 
results presented here about the proton and electron heating 
rates depend crucially on the use of the mean radial trends 
in the plasma parameters. Further work is needed to ensure 
that these results are valid for the totality of fast solar wind 
streams. Other approaches to estimating the effect of vari- 
ability of local parameters on the heating rates include those 
of Smith et al. (2006a) and Breech et al. (2008). 

It is useful to compare the empirically derived heating rates 
to those expected from a von Karman similarity analysis of 
the dissipation of an MHD turbulent cascade (e.g., Hossain et 
al. 1995). There are several pieces of evidence that suggest 
turbulence to be responsible for the in situ plasma heating 
(see, e.g., Coleman 1968). The most direct evidence is the 
fact that T p in the solar wind is positively correlated with the 
overall amplitude of the turbulent fluctuations (Grappin et al. 
1990; Vasquez et al. 2007). A more thorough discussion of 
the expected turbulent heating rate is given by, e.g., Breech et 
al. (2008, 2009). Roughly, though, we can estimate this rate 
as 

p(ZlZ. + Z 2 _Z + ) 

Ulmb ~ " 



A 



(15) 



where the mass density poc r and the transverse correlation 
length A^ oc r I//2 . Note that the heating rate Q tm b depends on 
the cross helicity of the fluctuations, which is defined by the 
ratio a c = (Zl-Z 2 )/(Zl + Z 2 ), and for highly Alfvenic states 
where a c — ► ±1, we see that gturb — > 0. 

The above phenomenological form for the turbulent heat- 
ing rate has been found to be consistent with numerical sim- 
ulations of strong MHD turbulence in a background magnetic 
field (e.g., Dobrowolny et al. 1980; Hossain et al. 1995; Zhou 
& Matthaeus 1990; Matthaeus et al. 1999a; Oughton et al. 
2001; Dmitruk et al. 2001, 2002). The success of transport 
theories in accounting for turbulence properties observed by 
Voyager, Ulysses, and Pioneer data (Zank et al. 1996; Smith 
et al. 2001; Breech et al. 2008) also provide empirical sup- 
port for the use of Qtmb in the form given by equation ([TBI . 



The radial scaling Aj_ oc r 1 ' 2 has been measured in the helio- 
sphere by, e.g., Bruno & Dobrowolny (1986) and Smith et al. 
(2001). The measured Elsasser amplitudes Z± exhibit a com- 
plex radial dependence that depends on various properties, 
but it is possible to approximate them reasonably well with 
WKB Alfven wave action conservation (Zank et al. 1996). 
Assuming spherical symmetry — which is appropriate for the 
high-speed wind at high latitudes — this radial dependence is 
Z± oc r" 1 / 2 . Putting these together we find that 2mrb oc r" 4 , 
which explains how the above approximations lead to the 
corresponding curves in Figure 2 being constant (see also 
Dmitruk et al. 2002; Cranmer & van Ballegooijen 2005). 

To normalize Qtm-b to the values shown in Figure 2, we ap- 
plied measured plasma properties at 1 AU, such as the proton 
number density (2.5 cm" 3 ) and the Elsasser amplitudes cor- 
responding to outwardly propagating (Z + ss 61 km s -1 ) and 
inwardly propagating (Z_ » 26 km s" 1 ) Alfven waves. The 
latter two quantities combine to give a mean Elsasser ampli- 
tude Z w 47 km s -1 at 1 AU and a normalized cross helicity 
<r c «s 0.7 (see, e.g., Bavassano et al. 2000; Dasso et al. 2005; 
Cranmer & van Ballegooijen 2005; Breech et al. 2008, 2009). 
The remaining free parameter in equation (IT~5b is the value 
of the transverse correlation length at 1 AU, which we vary 
between 0.01 and 0.2 AU. Figure 2 shows that the empirical 
heating rates for the fast wind are consistent with Aj_ w 0.02- 

0. 1 AU. Breech et al. (2008) found that this same range of 
correlation lengths applies to the fast wind at 1 AU as well. 6 

Figures 3 and 4 show the ratio of proton heating to the total, 

1. e., Q p /{Q P + Q e ), for a number of different calculations. This 
fraction is denoted f p by Breech et al. (2009). In all panels, 
a baseline model is shown for comparison that was computed 
with u = 700 km s" 1 and the standard choices for other param- 
eters as described above. Figure 3a shows the dependence of 
this ratio on the assumed value of the wind speed. The ratio 
at r < 1 AU appears to be insensitive to the wind speed be- 
cause both Q p and Q e vary linearly with u at these distances 
(see also Fig. 2). The increased spread in the ratio at larger 
distances is the result of the electron heat conduction having 
a larger relative impact on Q e as the wind speed changes. 

Figure 3b shows the dependence of the proton-to-total heat- 
ing ratio on varying several other assumptions of the standard 
model (all keeping u = 700 km s" 1 ). We computed trial ver- 
sions of the proton temperature fitting curve with the standard 
Ulysses measurements replaced by either the lower or upper 
limit measurements (r sma n or 7i aree as described in § 2). We 
also tried replacing the fitting curve for the electron heat con- 
duction flux (eq. 0) with the simple collisionless expression 
(eq. 0) and a constant coefficient a e = 1.05. This was seen 
in Figure lb to possibly be a reasonable description of the 
measured heat fluxes. 

It is noteworthy that the use of the collisionless model for 
qn e shows a rather extreme departure from the other results, 
with the ratio Q p /(Q p + Q e ) monotonically decreasing as a 
function of increasing distance. This difference arises because 
the collisionless heat flux remains steeper than r~ 2 in the in- 
ner heliosphere, and thus its divergence (see the last term in 
e Q- 0) gives a positive contribution to the electron heating at 
all distances. In a sense, the use of the collisionless heat flux 

6 Multispacecraft measurements (e.g., Matthaeus et al. 2005) show the 
measured correlation scale at 1 AU to be somewhat smaller than the values 
in Figure 2. This can be reconciled by using a va lue of order 0.1-0.2 for the 
von Karman proportionality constant in equation H5\ , instead of the value of 
1 assumed here (see also Breech et al. 2009). 
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FIG. 3. — Ratios of the proton heating rate to the total (proton + electron) 
heating rate. The standard model with wind speed u = 700 km s is shown in 
both panels (thick red line), (a) Variation of the wind speed between 600 and 
800 km s _1 in 50 km s _I increments (see labels), (b) Models computed with 
r sma u (dotted line) and 7] a „e (dashed line) for the SWOOPS proton temper- 
atures, and with equation (7) for qu e instead of the empirical fit (dot-dashed 
line). 

is consistent with the need for heat to be conducted outwards 
from heights below our lower boundary at 0.3 AU. The exis- 
tence of this additional conductive heating below 1 AU means 
that the deposited heating rate Q e need not be as large as it 
would be otherwise, and thus the protons end up contributing 
more to the total heating. On the other hand, when using the 
least-squares fit to the measured heat fluxes, the slope of qn e 
is flatter than r~ 2 in the inner heliosphere and conduction acts 
to cool the electrons at these distances. The derived value of 
Q e must then be larger, which gives the rough equipartition 
between protons and electrons. This is consistent with a more 
local deposition of heat to the electrons that conducts both up 
and down from some point around 0.75 AU. We should em- 
phasize, however, that the collisionless expression is based on 
relatively simple theoretical scalings and it is far from being 
a rigorous "prediction" for the electron heat flux in the solar 
wind. The discrepancy seen in Figure 3b suggests that the use 
of a constant a e coefficient in equation (JTJ should probably 
not be considered a robust approximation. 
Figure 4a shows the result of varying the Coulomb collision 
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FIG. 4. — Same as Figure 3, but with other parameters varied. (a)Aseriesof 
models with anomalously strong Coulomb collisions, with a range of constant 
multipliers / to the classical collision frequency (see labels), (b) Models 
computed over a range of colati tude s 9 in the heliosphere, which affects the 
Parker spiral angle <!> (see eq. 1101 ). The standard model with wind speed 
u = 700 km s" 1 is shown in both panels (thick red line). 



rate by multiplying v pe (eq. fT3l ) by a range of constant mul- 
tipliers /. As mentioned above, the curves for / = and / = 1 
would be indistinguishable and the former is not shown. As / 
is increased, collisions attempt to equalize T p and T e at a faster 
rate. The fraction of heat going to the protons must increase 
in order to maintain the specified T p > T e . Note, however, that 
above a certain value of / there would need to be net electron 
cooling, or essentially Q p > (Q p + Q e ), in order to maintain 
the measured temperature difference. This allows us to rule 
out values of / larger than about 200 over most heliocentric 
distances. Figure 4b shows how the relative heating changes 
as the assumed colatitude 9 (and thus the Parker spiral angle 
is varied. Models computed closer to the ecliptic plane ex- 
hibit a local maximum in Q p /(Q P + Q e ) at smaller heliocentric 
distances. 

Given the uncertainties in many of the input parameters to 
these models, the quantitative results regarding the relative 
partitioning of Q p and Q e are also somewhat uncertain. How- 
ever, there does appear to be a preponderance of evidence for 
the validity of two qualitative statements: (1) In the inner he- 
liosphere (r < 1 AU), there appears to be rough equipartition 
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between proton heating and electron heating. (2) As helio- 
centric distance is increased from 0.3 to 5 AU, the relative 
fraction of proton heating increases to noticeably dominate 
the total heating rate. 

4. EXPECTATIONS FROM LINEAR ALFVEN WAVE DAMPING 

It is illustrative to compare the empirically derived parti- 
tioning between proton and electron heating with theoretical 
predictions. As a starting point for future work in this direc- 
tion, we computed some extremely simple estimates for the 
wavenumber dependence of a turbulent power spectrum of 
Alfvenic fluctuations. These were then coupled to a general 
Alfven wave dispersion relation that allowed us to compute 
the contributions to proton and electron heating. Although it 
is well-known that strong MHD turbulence is far from "wave- 
like" (i.e., one might expect that a coherent wave survives 
for only about one period before nonlinear processes transfer 
its energy to smaller scales), there is a long history of using 
damped linear wave theory to study the small-scale dissipa- 
tion of a cascade (see, e.g., Eichler 1979; Quataert 1998; Lea- 
mon et al. 1999; Quataert & Gruzinov 1999; Marsch & Tu 
2001; Cranmer & van Ballegooijen 2003; Gary & Borovsky 
2004, 2008). A typical justification of this approach is that the 
amplitudes of magnetic fluctuations in the dissipation range 
tend to be extremely small ({SB 2 ) -C B^); see also Spangler 
(1991) and Miller etal. (1996). 

We assumed that the wavenumber dependence of the power 
spectrum of magnetic fluctuations Pg scales as kT 1 ! 2 , where 
k is the magnitude of the wavenumber. The power Pg(fc) is 
defined such that its integral over the full (three-dimensional) 
wavenumber space gives the total magnetic energy density of 
fluctuations. The exponent of -7/2 corresponds to an expo- 
nent of -3/2 for a corresponding one-dimensional isotropic 
spectrum (i.e., the latter measuring power in spherical shells 
in fc-space). This exponent has been proposed for various the- 
ories of MHD turbulence (e.g., Iroshnikov 1964; Kraichnan 
1965; Nakayama 2001). This value is also close to the ex- 
ponent of -10/3 that has been predicted for a purely perpen- 
dicular cascade having a one-dimensional power spectrum of 

k~f /3 (see Kolmogorov 1941; Fyfe et al. 1977; Higdon 1984; 
Goldreich & Sridhar 1995; Cho & Vishniac 2000; Boldyrev 
2005; Horbury et al. 2008). 

We also assumed that the above wavenumber dependence 
applies in restricted ranges of O^b, which is defined as the an- 
gle between the wavenumber vector and the background mag- 
netic field direction. Specifically, we constructed three mod- 
els with different angular distributions of wave power: (1) a 
"slab" spectrum, with nearly parallel-propagating waves fill- 
ing the region < O^b < 5°, (2) a "two-dimensional" (2D) 
spectrum, with nearly perpendicularly-propagating waves fill- 
ing the region 85° < O^b < 90°, and (3) an isotropic spectrum 
with all values of O^b having the same wavenumber depen- 
dence. 

The second case given above — i.e., a turbulent cascade 
mainly in the k± direction — uses a range of angles that is ac- 
tually quite broad in comparison to theoretical expectations 
of the so-called "critical balance" proposed by Goldreich & 
Sridhar (1995). This condition is defined by an equivalence 
between the time scales of Alfven wave propagation along the 
field and nonlinear energy transfer perpendicular to the field. 
According to Goldreich & Sridhar (1995), the majority of the 
power in strong MHD turbulence should have angles between 
the critical balance angle 8 cr [ t and 90°. At large perpendicular 



wavenumbers k± (i.e., at the onset of kinetic Alfven wave dis- 
persion and substantial Landau damping), we estimated that 
9 cr i t decreases slightly from approximately 89.8° at 0.3 AU to 
88.9° at 5 AU. These values were computed from the existing 
models of Cranmer & van Ballegooijen (2003, 2005) and the 
measured plasma parameters given in § 2. 

We computed the proton and electron heating rates for each 
kind of spectrum using the quasi-linear framework developed 
by Quataert (1998), Quataert & Gruzinov (1999), Marsch & 
Tu (2001), and Cranmer & van Ballegooijen (2003). The heat- 
ing rates are given by integrals over vector wavenumber k of 
the form 

q s = P J d\p ot (m^, (i6) 

where s = p,e denotes the particle type of interest and P tot w 
2P B is the total energy spectrum of fluctuations. The species- 
dependent amplitude damping rates w,^ were defined by Cran- 
mer & van Ballegooijen (2003) to be weighted by Lan- 
dau/cyclotron resonance functions that take account of the 
particle kinetic motions. For a proton-electron plasma, the 
sum LOi )P + u)i >e gives the absolute value of the total linear 
damping rate (i.e., the imaginary part of the wave frequency). 
The resonance functions depend on the dispersive properties 
of the Alfven waves, which were computed using the warm- 
plasma Vlasov-Maxwell dispersion code of Cranmer & van 
Ballegooijen (2003). These functions also depend on the 
shape of the particle velocity distributions, which were as- 
sumed here to be Maxwellian. 

In order to compute the relative heat dissipated by the pro- 
tons versus that dissipated by the electrons, we did not need 
to specify an absolute normalization for the power spectrum. 
However, in practice, it was useful to cut off the spectrum 
at low wavenumbers using an "outer scale" lower limit k out 
that scaled inversely with the turbulent correlation length dis- 
cussed above; i.e., k ou i oc oc r 2 . The upper limits for the 
parallel and perpendicular wavenumbers (k\\,k±) were found 
numerically by the Vlasov-Maxwell code as locations where 
the Alfven-wave solution branches ceased to exist (see Stix 
1992; Cranmer & van Ballegooijen 2003). 

The dispersive properties of Alfven waves depend strongly 
on the value of the plasma (3, which is defined here as the ratio 
of proton plasma pressure to the magnetic pressure (see Gary 
& Borovsky 2004, 2008). For the high-latitude heliocentric 
distances considered here, we used a very simple monopolar 
scaling for the magnetic field strength, 

BoW «2.5xlO" 5 ( T ^) 2 G (17) 

and (3 = %Ttn p kv,T p /B 2 y Using the empirically constrained den- 
sities and temperatures discussed above, the value of (3 in- 
creases monotonically from 1.1 (at 0.3 AU) to about 37 (at 
5 AU). Because calculating the dispersive properties of the 
waves on a fine two-dimensional grid in wavenumber is com- 
putationally expensive, we created only 7 grids having (3=1, 
1.8, 3.4, 6.4, 12, 22, and 40. We computed full runs of Q p 
and Q e versus heliocentric distance using each of the 7 grids, 
and then we interpolated between these results (using the em- 
pirical plasma (3 at each distance) in order to determine the 
appropriate solutions. 

Figure 5 shows the results of this procedure for the slab, 
2D, and isotropic spectra, and we compare the proton heat- 
ing fractions Q p /(Qp + Q e ) to the standard empirical curve 
from Figures 3 and 4. The curve denoting the slab spectrum 
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FIG. 5. — Comparison of the standard empirical proton-to-total heating 
ratio (solid line) to theoretical ratios computed from the linear damping of 
Alfven wave spectra: parallel-propagating "slab" modes only (dotted line), 
perpendicular "2D" modes only (dashed line), and an isotropic spectrum 
(dot-dashed line). 

is relatively simple to understand. In a collisionless plasma, 
parallel-propagating Alfven waves dissipate primarily by cy- 
clotron resonance with positive ions (see Hollweg & Isenberg 
2002). For the values of (3 used in the warm-plasma disper- 
sion code, the largest wave frequency for parallel-propagating 
waves was only about 0.3 times the local proton gyrofre- 
quency. Thus, since the model did not contain any other ions, 
only the protons felt the cyclotron resonance and were re- 
sponsible for almost all of the dissipation. There was some 
extremely weak Landau damping in the slab model that con- 
tributed about 0.1% of the total heating to the electrons. 

The curve denoting the 2D spectrum (i.e., fluctuations per- 
pendicular to the background magnetic field) always shows 
a roughly equal partitioning of energy between protons and 
electrons. However, the partition fraction oscillates slightly 
up and down several times between 0.3 and 5 AU. These os- 
cillations do not seem to be numerical artifacts, and their ori- 
gin is not entirely clear. For the range of plasma j3 found 
in these models, there is never a single physical process that 
dominates the collisionless damping of highly oblique Alfven 
waves. The mechanisms of electron Landau damping, proton 
Landau damping, and proton transit-time damping are all of 
comparable magnitude (e.g., Gary & Borovsky 2004, 2008). 
Thus, we refrain from attributing the mild variations around 
the mean equipartitioning (Q p « Q e ) to any specific physical 
process. We propose to study this complex interplay of colli- 
sionless dissipation mechanisms in future work. Finally, the 
curve denoting an isotropic fluctuation spectrum tends gives 
an intermediate amount of heat partitioning between the slab 
and 2D cases, as one would expect. 

By comparing the theoretical and empirical curves in Figure 
5, we can make some preliminary suggestions about how the 
anisotropy of the actual turbulent power spectrum (or at least 
its large-wavenumber dissipation range) varies between 0.3 
and 5 AU. An interpretation based on Figure 5 is that the fluc- 
tuations are close to 2D at the smallest heliocentric distances, 



and that they evolve to be more isotropic as distance is in- 
creased. This makes heuristic sense from at least one limited 
viewpoint: i.e., that as the plasma (3 increases and the effect 
of the magnetic field becomes less important to the plasma, 
the turbulence may want to behave in an increasingly hydro- 
dynamic (i.e., isotropic) fashion. Note, however, that this con- 
clusion is predicated on a wide range of assumptions about the 
models. Most importantly, it depends on the turbulent fluctu- 
ations being purely in the Alfven mode, and not having, say, 
fast-mode/whistler type dispersive properties. This conclu- 
sion also depends on the turbulent fluctuations damping like 
linear waves, despite our knowledge that nonlinear features 
such as current sheets may play an important role (Dmitruk et 
al. 2004). 

We should also note that if the fit to the electron heat con- 
duction (eq. 0) is incorrect, and instead the collisionless ex- 
pression of Hollweg ( 1 974, 1 976) is taken to be valid, then the 
above conclusions may be reversed. A comparison between 
the heating partition curve that was computed for a e = 1.05 
(from Fig. 3b) with the theoretical results in Figure 5 would 
imply that the turbulence evolves from a slab or isotropic 
spectrum at 0.3 AU to a more 2D spectrum at larger distances. 
This serves to emphasize how much our conclusions about 
proton and electron heating depend on the accurate knowl- 
edge of the radial dependence of qn e in the fast solar wind. 

There has been a great deal of work done to measure the 
anisotropy of MHD turbulence in the solar wind. Ultimately, 
the empirical data for both the particle heating and the fluctua- 
tion spectra should be combined to provide the tightest possi- 
ble constraints on theoretical models. At 1 AU, the dominant 
component of the turbulence does seem to be perpendicular 
or 2D (e.g., Bieber et al. 1996; MacBride et al. 2008). How- 
ever, when analyzing the radial dependence of the anisotropy, 
there have been differing answers to some of the key ques- 
tions. Some studies that compared the power between the 
minimum and maximum variance directions concluded that 
the turbulence becomes more isotropic as distance increases 
(Klein et al. 1991; Horbury et al. 1995). Other studies found 
that the perpendicular component becomes more dominant at 
larger distances (Bavassano et al. 1982; Neugebauer 2004). 
Such opposite conclusions may be the result of mixing to- 
gether regions that contain different amounts of fast and slow 
wind (see, e.g., Dasso et al. 2005) or using different ranges 
of frequency. Also, the variance anisotropy does not nec- 
essarily vary in the same way as the spectral (wavenumber) 
anisotropy. These results may also depend on subtleties of the 
analysis technique, such as how the minimum-variance direc- 
tion is interpreted (Smith et al. 2006b; Tessein et al. 2009) 
or how effects such as intermittency are treated (Bruno et al. 
1999). 

5. CONCLUSIONS AND FUTURE PROSPECTS 

The goal of this paper has been to compute empirical esti- 
mates for the rates of proton and electron heating in the fast 
solar wind between 0.29 and 5.4 AU. A key new aspect of this 
work has been the incorporation of the measured electron heat 
conduction flux, which dominates the internal energy balance 
of electrons in much of the heliosphere. We conclude that the 
protons receive about 60% of the total plasma heating in the 
inner heliosphere, and that this fraction increases to approx- 
imately 80% by the orbit of Jupiter. These results are con- 
firmed by the independent analyses of Leamon et al. (1999) 
and Breech et al. (2009), who found that a similar proton-to- 
total heating fraction of ^60% is consistent with theoretical 
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models of MHD turbulent dissipation. 

Various uncertainties in the measurements affected our re- 
sults. For example, the rates of proton and electron heating 
(Q p and Q e ) were calculated from mean radial trends in the 
data (as exemplified by the least-squares fits given in § 2), 
and they do not describe the substantial measured spread in 
the data about the mean values. Also, the uncertainty in 
Q e increases with heliocentric distance because the electron 
heat conduction flux becomes more important — in a relative 
sense — as distance increases. Small uncertainties in q\^ e or its 
radial slope thus have a larger impact on the derived heating 
rates at distances greater than ^2 AU. In any case, the general 
techniques presented in this paper are applicable to any future 
measurements that would improve on our current knowledge 
of the solar wind plasma properties. 

It is interesting to compare the interplanetary proton and 
electron heating rates with those inferred for the solar corona. 
The observational consensus since the late 1990s has been that 
coronal holes undergo preferential proton heating in a similar 
sense as high-speed wind streams in the heliosphere (see, e.g., 
Kohl et al. 1997, 1998, 2006; Wilhelm et al. 2007). Recently, 
Figure 6 of Landi & Cranmer (2009) summarized a range of 
measurements of T p and T e in polar coronal holes. Although 
one generally finds T p > T e , the measurements of T p and T e 
do not yet fully overlap with one another in heliocentric dis- 
tance. Improved measurements are needed in order to better 
constrain the proton and electron heating rates in the corona. 
However, given that the corona and the heliosphere are so dif- 
ferent in density, Alfven speed, and plasma (3, it is somewhat 
surprising that the proton-electron heat partitioning in these 
regions may be so similar in character. 

An improved understanding of the so-called two-fluid na- 
ture of solar wind plasma is an important ingredient in produc- 
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electron heating in global models may be a key to improving 
their physical realism and predictive accuracy. 
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